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Abstract 

The phase transition of hadronic to quark matter at high baryon and isospin 
density is analyzed. Nonlinear relativistic mean field models are used to describe 
hadronic matter, and the MIT bag model is adopted for quark matter. The bound- 
aries of the mixed phase and the related critical points for symmetric and asymmet- 
ric matter are obtained. Isospin effects appear to be rather significant. The binodal 
transition line of the (T, ps) diagram is lowered to a region accessible through heavy 
ion collisions in the energy range of the new planned facilities, e.g. the FAIR/NIC A 
projects. Some observable effects are suggested, in particular an Isospin Distillation 
mechanism. The dependence of the results on a suitable treatment of isospin con- 
tributions in effective QCD Lagrangian approaches, at the level of explicit isovector 
parts and/or quark condensates, is critically discussed. 
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1 Introduction 



Several suggestions are already present about the possibility of interesting 
isospin effects on the transition to a mixed hadron-quark phase at high baryon 
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Fig. 1. Zero temperature EoS of Symmetric/Neutron Matter: Hadron (NLp), solid 
lines, vs. Quark (MIT-Bag), dashed lines. an,Q represent the isospin asymmetry 
parameters respectively of the hadron, quark matter: ctu,Q = 0, Symmetric Matter; 
a H,Q = I, Neutron Matter. 

density [T1I21I3] . This seems to be a very appealing physics program for the 
new facilities, FAIR at GSI-Darmstadt [I] and NICA at JINR-Dubna [5], 
where heavy ion beams (even unstable, with large isospin asymmetry) will be 
available with good intensities in the 1-30 AGeV energy region. 

The weak point of those predictions is the lack of a reliable equation of state 
that can describe with the same confidence the two phases, hadronic and 
deconfined. On the other hand this also represents a strong theory motivation 
to work on more refined effective theories for a strong interacting matter. 

The main qualitative argument in favor of noticeable isospin effects on the 
hadron-quark transition at high density can be derived from the Figfl], where 
we compare typical Equations of State (EoS) for Hadron (Nucleon) and Quark 
Matter, at zero temperature, for symmetric (a = (p„ — p p )/pb = —PzIpb = 
0.0) and neutron matter (a = 1.0), where p n>p are the neutron/proton densities 
and Pb = Pn + pp the total baryon density. 

For the Hadron part we use a Relativistic Mean Field (RMF) EoS ([SHE]) 
with non-linear terms and effective p — meson coupling for the isovector part, 
largely used to study isospin effects in relativistic heavy ion collisions [3f8] . 
In order to keep a smooth flow of the physics points in the discussion, details 
about the adopted effective nucleon-meson Lagrangians are presented in the 
Appendix A. 
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The energy density and the pressure for the quark phase are given by the MIT 
Bag model [9] (two-flavor case) and read, respectively: 
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k 2 + m 2 (f q + f q ) + B , 
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d 3 k k 2 



=(f g + f q )-B, 
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(2tt) 3 ^k 2 + m\ 



where B denotes the bag constant (the bag pressure), taken as a rather stan- 
dard value from the hadron spectra [B = 85.7 MeV fm~ 3 , no density de- 
pendence), m q are the quark masses (m u = rrid = 5.5 MeV choice), and f q , 
f q represent the Fermi distribution functions for quarks and anti-quarks. The 
quark number density is given by 



The transition to the more repulsive quark matter will appear around the 
crossing points of the two EoS. We see that such crossing for symmetric matter 
(an — olq — 0.0) is located at rather high density, ps — 7po, while for pure 
neutron matter (an = olq = 1.0) it is moving down to about three times po- 
Of course the Fig.l represents just a simple energetic argument to support the 
hadron-quark transition to occur at lower baryon densities for more isospin 
asymmetric matter. In the rest of the paper we will rigourously consider the 
case of a first order phase transition in the Gibbs frame for a system with 
two conserved charges (baryon and isospin), in order to derive more detailed 
results. Since the first order phase transition presents a jump in the energy, 
we can expect the mixed phase to start at densities even before the crossing 
points of the Fig.l. The lower boundary then can be predicted at relatively low 
baryon densities for asymmetric matter, likely reached in relativistic heavy ion 
collisions. Moreover this point is certainly of interest for the structure of the 
crust and the inner core of neutron stars, e.g. see refs. [TOfTl] and the review 
\1'2\ . The ref. [TD] is particularly interesting since similar results are obtained 
with rather different hadronic approaches, the RMF and the non-relativistic 
Brueckner-Hartree-Fock (BHF) theory. 

We finally note that the above conclusions are rather independent on the 
isoscalar part of the used Hadron EoS at high density, that is chosen to be 
rather soft in agreement with collective flow and kaon production data [T3"fT3] . 
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In the used Bag Model no gluon interactions, the a s -strong coupling parame- 
ter, are included. We remark that this in fact would enhance the above effect, 
since it represents an attractive correction for a fixed B-constant, see |15j . A 
reduction of the Bag-constant with increasing baryon density, as suggested by 
various models, see ref. [TO], will also go in the direction of an "earlier" (lower 
density) transition, as already seen in ref. [2]. At variance, the presence of ex- 
plicit isovector contributions in the quark phase could play an important role, 
as shown in the following also for other isospin properties inside the mixed 
phase. 



2 Isospin effects on the Mixed Phase 



We can study in detail the isospin dependence of the transition densities [Tf2"f5] . 
The structure of the mixed phase is obtained by imposing the Gibbs conditions 
[T6] for chemical potentials and pressure and by requiring the conservation of 
the total baryon and isospin densities: 



M,P%,T) = » Q b{p% PIT), 

pnpB,P^T)=^(p Q B , P iT) , 

P H {T){p H B) p^ T) = pQ(T)(p Q B ,plT), 
Pb = (1- x)Pb + XP% > 

Ps = (1 - x)Pz + XP% , (4) 



where \ is the fraction of quark matter in the mixed phase and T is the 
temperature. 

The consistent definitions for the densities and chemical potentials in the two 
phases are given by : 



PB = P P + Pn, P3 = Pp~ Pn , 

H Pp Pn H Pp Pn (r\ 

Pb — 2 ' ^ 3 — 2 ' 
for the Hadron Phase and 



Pb = 7;(pu + pd) 



,Q 



/¥ = Pn - Pd , 
Pu — Pd 



(6) 
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for the Quark Phase. 



The related asymmetry parameters are: 



a 




Pn ~ Pp 
Pn + Pp 



a 



Q 




3 



Pd ~ Pu 
Pd + Pu 



(7) 



Nucleon and quark chemical potentials, as well as the pressures in the two 
phases, are directly derived from the respective EoS. 

In this way we get the binodal surface which gives the phase coexistence region 
in the (T,p B , p 3 ) space. For a fixed value of the total asymmetry ax = —pz/ps 
we will study the boundaries of the mixed phase region in the (T,ps) plane. 
Since in general the charge chemical potential is related to the symmetry term 
of the EoS, [8], ^3 = 2E sym (pB) — , we expect critical and transition densities 
rather sensitive to the isovector channel in the two phases. 

In the hadron sector we will use the Non-Linear Relativistic Mean Field mod- 
els, [7,8,3J, with different structure of the isovector part, already tested to 
describe the isospin dependence of collective flows and meson production for 
heavy ion collisions at intermediate energies, [17,18,19]. We will refer to these 
different Iso-Lagrangians as: i) NL, where no isovector meson is included and 
the symmetry term is only given by the kinetic Fermi contribution, ii) NLp 
when the interaction contribution of an isovector-vector meson is considered 
and finally iii) NLpS where also the contribution of an isovector-scalar meson 
is accounted for. See details in Appendix A and refs. [7118113] . 

We will look at the effect on the hadron-quark transition of the different stiff- 
ness of the symmetry term at high baryon densities in the different parametriza- 
tions. As clearly shown in Appendix A, where a rather transparent form for 
the density dependence of the symmetry energy in RMF approaches is dis- 
cussed, the potential part of the symmetry term will be proportional to the 
baryon density in the NLp choice and even stiffer in the NLpS case. 

As already mentioned, in the quark phase we use the MIT-Bag Model, where 
the symmetry term is only given by the Fermi contribution. The Bag param- 
eter B is fixed for each baryon density to a constant, rather standard, value 
B l l A = 160MeV, corresponding to a Bag Pressure of 85.7 MeV fm~ 3 . 

In general for each effective interactive Lagrangian we can simulate the solu- 
tion of the highly non-linear system of Eqs.flU), via an iterative minimization 
procedure, in order to determine the binodal boundaries and the Critical End 
Point (CEP) (T c , pf ) of the mixed phase. 

A relatively simple calculation can be performed at zero temperature. The 
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Fig. 2. Dependence on the Hadron Symmetry Energy of the Lower (x = 0.0) 
and Upper {x = 1.0) Boundaries of the Mixed Phase, at zero temperature, vs. the 
asymmetry parameter. Quark EoS: MIT bag model with i? 1//4 =160 MeV . 



isospin effect (asymmetry dependence) on the Lower (x = 0.0) and Upper 
(x = 1.0) transition densities of the Mixed Phase are shown in Figj2] for 
various choices of the Hadron EoS. The effect of a larger repulsion of the 
symmetry energy in the hadron sector, from NL to NLp and to NLpS, is 
clearly evident on the lower boundary with a sharp decrease of the transition 
density even at relatively low asymmetries. 

Typical results for isospin effects on the whole binodal "surface" are presented 
in Figj3]for symmetric and asymmetric matter. For the hadron part we have 
started from a NLp effective Lagrangian very close to other widely used rela- 
tivistic effective models, e.g. see the GM3 of ref.J20] and the NL3 interaction 
of P.Ring and collaborators [21], which has also given good nuclear structure 
results, even for exotic nuclei. 

In the symmetric matter case the mixed phase is evaluated from the simpler 
Maxwell conditions. The results are shown in FigJUfor the same hadron and 
quark EoS's as in Fig 13] at temperatures T=0, 50 and 80 MeV. The equal 
chemical potential densities (intersection of the dotted line in the lower panel) 
must correspond to the equal pressure densities of the upper panels. We nicely 
see: i) at T=0 MeV the mixed phase is centered around pj p — 7.0 , exactly 
the a = crossing point of Fig]]], confirming our energetic argument about 
the transition location; ii) the size of the mixed phase is shrinking with tem- 
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Fig. 3. Binodal surface for symmetric (a = 0.0) and asymmetric (a = 0.2) matter. 
Hadron EoS from NLp interaction. Quark EoS: MIT bag model with 5 1 / 4 =160 
MeV. The grey region corresponds to the conditions reached in Heavy Ion Collisions 
simulations at few AGeV beam energies, see Section 3. 



T=0 MeV T=50 MeV T=80 MeV 




Fig. 4. Maxwell construction for symmetric (a = 0.0) matter at temperatures T=0, 
50 and 80 MeV. For Pressure and Chemical Potentials we use: Hadron EoS (black 
curves) from NLp interaction; Quark EoS (grey curves) from MIT bag model with 
S 1 / 4 =160 MeV. 
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perature, it is very narrow at T=50 MeV and finally at T=80 MeV we cannot 
have anymore a first order transition. In fact the Critical End Point is found 
at T c ~ hSMeV, p c /p ~ 3.8, P c ~ l20MeV/fm 3 and fi e ~ U30MeV, see 
also FigJHJ The result is dependent on the choice of the Bag constant, with 
increasing of the critical temperature with the Bag value due to the reduction 
of the pressure in the quark phase, while the chemical potentials are not af- 
fected. This point will not be further discussed since here the main focus is on 
the isospin dependence of the mixed phase at lower temperatures, that can be 
probed in heavy ion collisions at intermediate energies. 

As expected, the lower boundary of the mixed phase is mostly affected by 
isospin effects. In spite of the relatively small total asymmetry, a = 0.2, we 
clearly observe in FigJH] a shift to the left of the first transition boundary, 
in particular at low temperature, and an indication of a relatively "earlier" 
Critical End Point, around (T c = 50 - 55MeV,p c /p = 2.5 - 2.8). 

Actually we see from FigJ3]that for asymmetric matter, we are not able to ex- 
plicitly reach the CEP of the binodal surface. Since this feature appears also 
in the following figures we must add a general comment. From the present 
results we cannot exclude a CEP at higher temperature and smaller baryon 
density. However we have also to mention some numerical problems. The solu- 
tion of the Gibbs conditions, the highly non-linear system of Eqs.fll]), is found 
through an iterative multi-parameter minimization procedure (the Newton- 
Raphson method). When we are close to the Critical End Point, for each 
^-concentration the baryon and isospin densities of the hadron and quark 
phases become very similar and we start to have problems in finding a definite 
minimum. We are working on improving the numerical accuracy but in any 
case the trends towards the CEP are always quite clearly observed, as a severe 
reduction of the mixed phase. 

In the following we will concentrate on properties of the mixed phase mostly 
located at high density and relatively low temperature, well described in the 
calculation and within the reach of heavy ion collisions in the few AGeV range, 
see Section 3. 



2.1 Inside the Mixed Phase of Asymmetric Matter 

For a = 0.2 asymmetric matter, in the FigsJSl we show also the (T,ps) 
curves inside the Mixed Phase corresponding to a 20% and 50% presence 
of the quark component (x = 0.2,0.5), evaluated respectively with the two 
choices, NLp and NLp5, of the symmetry interaction in the hadron sector. 
We note, as also expected from Figf2j that in the more repulsive NLp5 case 
the lower boundary is much shifted to the left. However this effect is not so 



8 




Fig. 5. Asymmetric a = 0.2 matter. Binodal surface and (T,ps) curves for various 
quark concentrations (x = 0.2, 0.5) in the mixed phase. Quark EoS: MIT Bag 
model with 5 1//4 =160 MeV. Hadron EoS: NLp Effective Interaction. 




Fig. 6. As in FigHJ for the NLpS Effective Interaction in the Hadron sector. 

evident for the curve corresponding to a 20% quark concentration, and almost 
absent for the 50% case. The conclusion seems to be that for a stiffer symmetry 
term in a heavy-ion collision at intermediate energies during the compression 
stage we can have more chance to probe the mixed phase, although in a region 
with small weight of the quark component. 
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Fig. 7. Asymmetric a = 0.2 matter. Binodal surface and (T, p^ , p^) curves for 
various quark concentrations in the mixed phase. Quark EoS: MIT Bag model 
with B 1 / 4 =160 MeV. Hadron EoS: NLp Effective Interaction. 
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Fig. 8. As in FigJTJ for the NLp8 Effective Interaction in the Hadron sector. 

In fact from the solution of the system Eq.Q we get the baryon densities 
PbiPb m the two phases for any x value. In the Figs. [7J [8] we present the 
results for the same weights 20%, 50% of the quark phase of the previous 
figures. The quark phase appears always with larger baryon density, even for 
the lowest value of the concentration. 

Can we expect some signatures related to the subsequent hadronization in the 
following expansion? 

An interesting possibility is coming from the study of the asymmetry a® in the 
quark phase. In fact since the symmetry energy is rather different in the two 
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Fig. 9. Quark asymmetry in the mixed phase vs. the quark concentration for asym- 
metric matter with T = and a = 0.2. NLp and NLp5 Effective Hadron Interac- 
tions are considered. Quark EoS: MIT bag model with Z? 1//4 =160 MeV. 

phases we can expect an Isospin Distillation (or Fractionation), very similar 
to the one observed in the Liquid-Gas transition in dilute nuclear matter 
[22ll23|8] . this time with the larger isospin content in the higher density quark 
phase. 

In Figj9] we show the asymmetry a® in the quark phase as a function of the 
quark concentration x f° r the case with global asymmetry a = 0.2 (zero tem- 
perature). The calculation is performed with the two choices of the symmetry 
term in the hadron sector. We see an impressive increase of the quark asymme- 
try when we approach the lower boundary of the mixed phase, even to values 
larger than one, likely just for numerical accuracy [21]. Of course the quark 
asymmetry recovers the global value 0.2 at the upper boundary \ = 1. A 
simple algebraic calculation allows to evaluate the corresponding asymmetries 
of the hadron phase. In fact from the charge conservation we have that for 
any x-mixture the global asymmetry a is given by: 

_ Ps (1 - x)a H , xa Q , oN 
a - 



PB (i-x) + x§ (i-x)4+x 



For any Xi from the calculated a Q of Figj9] and the Pb,Pb °f Fig s E El we 
can get the correspondent asymmetry of the hadron phase a H . For a 20% 
quark concentration we have an a®/a H ratio around 5 for NLp and around 
20 for NLp5, more repulsive in the isovector channel. It is also interesting to 
compare the isospin content N/Z of the high density region expected from 
transport simulations without the Hadron-Quark transition and the effective 
N/Z of the quark phase in a 20% concentration. In the case of Au + Au 
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(initial N/Z = 1.5) central collisions at lAGeV in pure hadronic simulations 
we get in the high density phase a reduced N/Z ~ 1.2 — 1.25 (respectively 
with NLp5 — NLp interactions) due to the fast neutron emission [r8"][T9"] . The 
corresponding isospin content of the quark phases is much larger, N/Z = 3 for 
NLp and N/Z = 5.7 for NLpS. This is the neutron trapping effect discussed 
in the Section 3. We could expect a signal of such large asymmetries, coupled 
to a larger baryon density in the quark phase, in the subsequent hadronization. 

We finally remark that at higher temperature and smaller baryon chemical 
potential (ultrarelativistic collisions) the isospin effects discussed here are ex- 
pected to vanish [25], even if other physics can enter the game and charge 
asymmetry effects are predicted also at p B = and T ~ 170 MeV [26|27] . 



Isospin in Effective Quark Models 

All the above results will be very sensitive to the explicit inclusion of isovec- 
tor interactions in effective non-perturbative QCD models in regions at high 
baryon chemical potentials. Unfortunately few attempts have been worked out 
for two main reasons: i) the difficulties of Lattice-QCD calculations at high 
baryon densities; ii) the main interest on the QGP phase transition at high 
temperature and small baryon chemical potentials, as probed in the expand- 
ing fireball of ultrarelativistic heavy ion collisions. A first approach can be 
supplied by a two-flavor Nambu-Jona Lasinio (NJL) model [28], which in 
fact describes the chiral restoration but not the deconfinement dynamics. The 
isospin asymmetry can be included in a flavor- mixing picture [2"9T3U] . corre- 
sponding to different couplings to the (u,d) quark-antiquark condensates. As a 
consequence we can have now a dependence of the constituent mass of a given 
flavor to both quark condensates. We devote the Appendix B to a detailed 
study of this isospin effects in the NJL chiral dynamics. 

Due to the scalar nature of the interacting part of the corresponding La- 
grangians only the quark effective mass dynamics will be affected. In the "re- 
alistic" small mixing case, see also [29.31J, we get a definite M* > splitting 
at high baryon density (before the chiral restoration). 

All that can indicate a more fundamental confirmation of the m* > m* split- 
ting in the hadron phase, as suggested by the effective QHD model with the 
isovector scalar 5 coupling, see [7]I5] . However such isospin mixing effect results 
in a very small variation of the symmetry energy in the quark phase, related 
only to the Fermi kinetic contribution. Moreover we remind that confinement 
is still missing in this NJL mean field approach. In any case there are exten- 
sive suggestions about a favored chiral symmetry restoration in systems with 
large neutron excess [55] . 
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More generally starting from the QCD Lagrangian one can arrive to an effec- 
tive color current-current interaction where an expansion in various compo- 
nents can provide isovector contributions, |34J. 

In this respect we remark another interesting "indirect" isospin effect, i.e. not 
directly related to isovector terms in the effective Lagrangian, related to the 
presence of quark condensates due to the attractive gluon interaction. We note 
that just a color-pairing effect in the two- flavor system (the 2SC phase [35] ) 
would imply a stiffer symmetry energy in the quark EoS since we have a larger 
attraction when the densities of up and down quarks are equal. A first study 
of the high density hadron-quark transition including such gluon correlation 
in the Bag model has been presented very recently in |36j . Now the symmetry 
energy difference between hadron and quark phases is partially reduced, at 
least at low temperatures, and consequently also the isospin effects discussed 
in detail in this work will be less robust, although still present. An interesting 
point is that in any case the quark phase is more bound due to the attractive 
gluon contribution. Hence the transition to the mixed phase will still appear at 
relatively low baryon densities, now for an "isoscalar" mechanism, within the 
reach of "low energy" heavy ion collisions, i.e. in the range of few AGeV. As 
an intuitive picture we can refer again to the Fig.l. Essentially the difference 
between the aQ = 0.0 and o.q — 1 curves is increasing but meanwhile both 
are decreasing. 

With increasing temperature the color pairing effect will be in general reduced, 
as confirmed in [37] in an extended NJL calculation, and so isospin effects, as 
discussed before, will be more relevant. All that, as well as the location of the 
Critical-End- Point, is naturally related to the used value of the superconduct- 
ing gap, opening new stimulating perspectives. In this sense new experiments 
on mixed phase properties observed with isospin asymmetric heavy ion colli- 
sions, as suggested in the final section, will be extremely important. 



3 Perspectives and suggested observables 

Based on the qualitative argument of the Introduction and on more detailed 
calculations in a first order phase transition scheme, we have predicted rather 
"robust" isospin effects on the hadron-quark transition at high baryon densi- 
ties, not depending on details of the EoS parametrizations in the hadron and 
quark phases. 

Our results seem to indicate a specific region where the onset of the mixed 
phase should be mainly located: 2 < Pb/po < 4, T < 50 — 60 MeV, for realistic 
asymmetries a ~ 0.2 — 0.3. A key question is if such a region of the phase space 
can be explored by means of Heavy- Ion-Collisions. In refs. [2*13] it is shown that 
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even collisions of stable nuclei at intermediate energies (E/A ~ 1 — 2 GeV) 
make available the pertinent (T,pQ,a) region where the phase transition is 
expected to occur. 

In this respect we can refer to the reaction 238 JJ + 238 JJ (average isospin asym- 
metry a = 0.22) at 1 AGeV that has been investigated in ref.p], using a 
consistent Relativistic Mean Field approach with the same interactions, for a 
semicentral impact parameter 6 = 7 fm, chosen just to increase the neutron 
excess in the interacting region. The evolution of momentum distribution and 
baryon/isospin densities in a space cell located in the cm. of the system has 
been also studied. After about 10 fm/c a local equilibration is achieved still 
in the compressed phase, before the fast expansion. We have a unique Fermi 
distribution and from a simple fit the "local" temperature can be evaluated. 
A rather exotic nuclear matter is formed in a transient time of the order of 
10 — 20 fm/c, with baryon density around 3 — 4p , temperature 50 — 60 MeV, 
and isospin asymmetry between 0.2 and 0.3, likely inside the estimated mixed 
phase region. We note that high intensity 238 U beams in this energy range 
will be available in the first stage of the FAIR facility [4.38] and also at JINR- 
Dubna in the Nuclotron first step of the NIC A project [25] . 

Which are the observable effects to look at if we enter and/or cross the mixed 
phase? 

As already stressed, a first expectation will be the Isospin Distillation effect, a 
kind of neutron trapping in the quark phase, supported by statistical fluctua- 
tions [2] as well as by a symmetry energy difference in the two phases, as dis- 
cussed in Section 2.1. In fact while in the pure hadron matter (neutron-rich) at 
high density we have a large neutron potential repulsion (in both NLp, NLp5 
cases), in the quark phase the rf-quarks see a much smaller symmetry repul- 
sion only due to the kinetic contribution from the Fermi gas. As a consequence 
while in a pure hadronic phase neutrons are quickly emitted or "transformed" 
in protons by inelastic collisions [19], when the mixed phase starts forming, 
neutrons are kept in the interacting system, in the quark phase which is also 
at larger baryon density, up to the subsequent hadronization in the expansion 
stage [3J. Observables related to such neutron "trapping" could be an inver- 
sion in the trend of the formation of neutron rich clusters, an enhancement 
of the production of isospin-rich nucleon resonances and subsequent decays, 
i.e. an increase of 7r~/7r + , K°/K + yield ratios for reaction products coming 
from high density regions, that could be selected looking at large transverse 
momenta, corresponding to a large radial flow. 

If such kinetic selection of particles from the mixed phase can really be suc- 
cessful also other potential signatures would become available. One is related 
to the general softening of the matter, due to the contribution of more degrees 
of freedom, that should show up in the damping of collective flows [39] . 
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The azimuthal distributions (elliptic flows) will be particularly affected since 
particles mostly retain their high transverse momenta escaping along direc- 
tions orthogonal to the reaction plane without suffering much rescattering 
processes. Thus a further signature could be the observation, for the selected 
particles, of the onset of a quark-number scaling of the elliptic flow: a property 
of hadronization by quark coalescence that has been predicted and observed 
at RHIC energies, i.e. for the transition at fiB = [ID]. 

We note that all the above results, on the Binodal Boundaries of the mixed 
phase and on the Isospin Distillation are sensitive to the symmetry term in 
the hadron sector, although the main isospin effects are present for all the 
parametrizations of the isovector interaction. At variance, for the quark sector 
the lack of explicit isovector terms could strongly affect the location of the 
phase transition in asymmetric matter and the related expected observables. 

In conclusion the aim of this work is twofold: 

• To stimulate new experiments on isospin effects in heavy ion collisions at 
intermediate energies (in a few AGeV range) with attention to the isospin 
content of produced particles and to elliptic flow properties, in particular 
for high-p t selections. 

• To stimulate more refined models of effective Lagrangians for non-perturbative 
QCD, where isovector channels are consistently accounted for and/or gluon 
correlations, leading to diquark condensates, can induce symmetry energy 
effects. 
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A Equation of state for hadronic matter 

Nonlinear (NL) relativistic mean field model with constant couplings 

A Lagrangian density with "minimal" meson channels and non-linear terms 
is used. The nuclear interaction is mediated by two isoscalar, the scalar a 
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and the vector u, and two isovector, the scalar S and the vector p, mesons. 
Non linear terms are considered only for the a contribution to account for 
the correct compressibility around saturation. Constant nucleon-meson cou- 
plings are used, chosen to reproduce the saturation properties and to repre- 
sent a reasonable average of the density dependence predicted by Relativistic 
Brueckner-Hartree-Fock calculations f£Tfl2] . see details in refs. [7|8] . 



L = #7/^ - (M - g a o - g 5 f ■ 8) - g^co" - g p ^f ■ b^tp 
+\{d^a - m 2 a a 2 ) - U(<j>) + \mlu^ + ±m% ■ 

+±{dj- d»5-m 2 P) - \F, V F^ - \g» v &\ (A.l) 



where F^ v = d^co v — d u u p , G pu = d^b u — d u b p , and the U(a) is the nonlinear 
potential of a meson : U(a) = |acx 3 + \ba^. 

The EoS for nuclear matter at finite temperature in the mean-field approxi- 
mation (MFA) is given by 



2E 



1 

2 m 



i=n,p 
2 



S^Etmm) + fi(k)) + \™W + U{a) 



+-— Pb + o^Ps + 



2mf s3 



(A.2) 



2m* PB + 2mf 3+ 2mf s3 > 



where E* = Jk 2 + M* 2 . The nucleon effective masses are defined as 

M* = M — g a (j) =f gs5 3 (— proton, + neutron). (A. 4) 

The field equations in the relativistic mean field (RMF) approach are 
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Oil Q bo Q (J / \ / A \ 

O = J(T j(J H ^{Psp + Psnj, (A.5) 



m„ m% m" 
w = -^P , (A.6) 

&o = ^Ps , (A.7) 

#3 = ^oiPsp - Psn) , (A.8) 

mi 



with the baryon density p = p^ = p p + p n and p% = p p — p n , p sp and p sn are 
the scalar densities for proton and neutron, respectively. The fi(k) and fi(k) 
in Eqs. (1A.2l (1A.3|) are the fermion and antifermion distribution functions for 



protons and neutrons {i = p,n): 



1 

m = i + exp{TO +/1 j)/r}- (A - 10) 



where the effective chemical potential p* is determined by the nucleon density 

Pi 

Pi = 2f ^(/*(*) -/*(*)) , (A.ll) 



while the scalar density p Sji , which gives the coupling to the scalar fields is 
given by 



r H 3 A- M* 

Ps, = 2j j^y 3 ^(Mk) + fm, (A.12) 

note the M* / E* quenching factor at high baryon density. 

The p* are related to the chemical potentials pi = de/dpi in terms of the 
vector meson mean fields by the equation 



Q 2 q 2 

Pp = P p + —PB + ^P3 

g 2 9 2 

Pn=P* n + ~^PB ~ (A.13) 



P 
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The baryon and isospin chemical potentials in the hadron phase can be ex- 
pressed in terms of the (p,n) ones as 



h _ Py j~ /A 



H 



2 



(A.14) 



In presence of the coupling to the two isovector p, 5-meson fields, the expression 
for the symmetry energy has a simple transparent form, see [71H3|8] : 



Now we easily see that in the NLpd choice we have a large increase of the 
symmetry energy at high baryon densities. The potential simmetry term is 
given by the combination [f p — fs(M*/E F ) 2 ] of the repulsive vector p and 
attractive scalar 5 isovector couplings. Thus, when the 6 is included we have 
to increase the p-meso coupling in order to reproduce the same asymmetry 
parameter a± at saturation. The net effect will be a stiffer symmetry energy at 
higher baryon densities due to the M*/E F quenching of the attractive part. 



Parameter determination 



The coupling constants are fixed from good saturation properties and 
from Dirac-Brueckner estimations, see the detailed discussions in refs. [7E31I8] . 
The isoscalar part of the EoS is chosen to be rather soft at high densities, see 
|44j . in order to satisfy the experimental constraints from collective flows and 
kaon production in intermediate energy heavy ion collisions [T3fT4"] . 

The coupling constants, f; L = gf/m 2 , i = a,u,p,5, and the two parameters of 
the a self-interacting terms : A = afg\ and B = b/g^ are reported in Table 
1. The a mass is fixed at 550MeV. The corresponding properties of nuclear 
matter are listed in Table 2. Here the binding energy is defined E/A — e/ p—M. 



7—i / \ 1 k fP 1 r p , \0~\ 

E sym {p) = 7^7 + M~ /*( pT) \P 



(A.15) 




where M* = M - g a a and Ep* = \jk 2 F + M* 2 . 



Table 1. Parameter set. 
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Parameter Set 


NLp 


NLp5 


U (fm 2 ) 


10.32924 


10.32924 


L (fm 2 ) 


5.42341 


5.42341 


f P (fm 2 ) 


0.94999 


3.1500 


fs (fm 2 ) 


0.000 


2.500 


A (fm- 1 ) 


0.03302 


0.03302 


B 


-0.00483 


-0.00483 



Table 2. Saturation properties of nuclear matter. 



Po (fm 3 ) 


0.16 


E/A (MeV) 


-16.0 


K (MeV) 


240.0 


E sym (MeV) 


31.3 


M*/M 


0.75 



We finally note that these Lagrangians have been already used for flow [T7] . 
pion production [18], isospin tracer (45] and kaon production [19] calculations 
for relativistic heavy ion collisions with an overall good agreement to data. 



B Nambu-Iona Lasinio model for asymmetric matter 

From the above discussion it appears (extremely) important to include the 
Isospin degree of freedom in any effective QCD dynamics. A first approach 
can be supplied by a two-flavor Nambu-Jona Lasinio model where the isospin 
asymmetry can be included in a flavor- mixing picture [29|30] . The lagrangian 
is given by 

L = L + L 1 + L 2 , (B.l) 
with Lq the free part 

L Q = ip(i $ — m)ip, 
and the two different interaction part given by 
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L 2 = G 2 {(#) 2 - (tfjrtfj) 2 - $i lh iP) 2 + (^t l5 r^) 2 }. (B.2) 

In the mean field approximation the new Gap Equations are Mj = — 
4G 1 $ i — 4G 2 $j, i 7^ j, (u, d) , where the $ Ui d =< m« >, < dd > are the two 
(negative) condensates which are given by 

*/ = -2iVc/ ^ ^{l - /-(T.ai/) - f + (T,Vf)}. (B.3) 
and m Mj d = m the (equal) current masses. 

Introducing explicitily a flavor mixing, i.e. the dependence of the constituent 
mass of a given flavor to both condensate, via G\ — (1 — p)(jQ) G% — @G@ we 
have the coupled equations 

M u = m- 4G $ U + 4/3G ($ u - $ d ), 

M d = m- AGoK + 4(1 - {3)G ($ U - $ d ). (B.4) 

For P = 1/2 we have back the usual NJL (M u = Md), while small/large mixing 
is for 0//3 ^ 1 respectively. The value of /3 has a consequence on the 
structure of the phase diagram in the region of low temperatures and high 
chemical potential. In fact as shown in [29|30] for /3 = there are two distinct 
phase transitions for the up quarks and for the down quarks, but for this value 
the interaction is symmetric under Ua(X) transformations and it is unrealistic. 
While for /3 > 0.1 the C/a(1) symmetry becomes explicitly broken and there 
is only a single first order phase transition. Realistic estimations of /3 fitting 
the physical //-meson mass give a value of /3 ~ 0.11 [29f3T] . 

In neutron rich matter | $ d | decreases more rapidly due to the larger p^, and 
so ($ n — <3>d) < 0. In the "realistic" small mixing case we will get a definite 
M u > Ma splitting at high baryon density (before the chiral restoration). This 
expectation is confirmed by a full calculation of the coupled gap equations 
with standard parameters [3J32] . All that can indicate a more fundamental 
confirmation of the m* > m* splitting in the hadron phase, as suggested by 
the effective QHD model with the isovector scalar 5 coupling, see [7f5] . 

However such isospin mixing effect results in a very small variation of the 
symmetry energy in the quark phase, still related only to the Fermi kinetic 
contribution. In fact this represents just a very first step towards a more 
complete treatment of isovector contributions in effective quark models, of 
large interest for the discussion of the phase transition at high densities. 
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